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Abstract 

We propose a model to represent the motility of social elements. The model is completely 
deterministic, possesses a small number of parameters, and exhibits a series of properties that are 
reminiscent of the behavior of comunities in social-ecological competition; these are: (i) similar 
individuals attract each other; (ii) individuals can form stable groups; (iii) a group of similar 
individuals breaks into subgroups if it reaches a critical size; (iv) interaction between groups can 
modify the distribution of the elements as a result of fusion, fission, or pursuit; (v) individuals can 
change their internal state by interaction with their neighbors. The simplicity of the model and 
its richness of emergent behaviors, such as, for example, pursuit between groups, make it a useful 
toy model to explore a diversity of situations by changing the rule by which the internal state of 
individuals is modified by the interactions with the environment. 
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I. INTRODUCTION 



Many living systems display collective organized displacements. Bird flocks, flsh schools, 
animal herds, bacterial colonies, and swarms are examples of such behavior. The interactions 
between the traveling individuals determine the collective motion. This kind of phenomenon 
been modeled by pb..e.t. Q. fl fl B Bfl. In the above .ode. all eleven, in 
the system are identical, but Shibata and Kaneko 9|] considered a coupled map gas in which 
each element is characterized by an internal state that evolves as a result of interactions with 
its neighbors and additionally, the forces between element pairs depend on their internal 
states. In [9| (see also [10|) the total force acting on a given element i is the sum of pair 
forces Fi j exerted by each of the elements j in its vicinity, and Fi j = —Fj i. This kind 
of reciprocity is appropriate to describe many physical systems, but in social systems the 
reaction of an individual to its environment is usually the result of a global evaluation instead 
of an evaluation of the sum of neighbors' actions. In other words, the reaction of a social 
individual i to the presence of individual j does not have the same intensity as and the 
opposite direction to the reaction of j to the presence of i. In particular, persecution and 
flight are not possible if Fij = —Fj^i. This kind of behavior requires a global evaluation of 
the environment before a decision to initiate a persecution or an escape. 

This paper is organized as follows. In Sec. II we present a simple model of interacting 
motile elements which display pattern dynamics that evoke some basic behaviors of social 
comunities. The interactions of the motile element have a flnite range and modify both 
their internal states and their positions. The functional form of the coupling of an element 
with its vicinity allows the formation of groups, fusion and fission of grou ps, and persecution 
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14, hj. In Sec. 



and flight of groups, behavior that is observed in many species 
III we analyze the properties of an isolated group as a function of its size and obtain the 
maximum density of elements for a compact conflguration of isolated groups. In Sec. IV we 
show detailed simulations in a two-dimensional space for three cases of internal dynamics 
corresponding to (i) a homogeneous and stationary internal state, (ii) two stationary but 
opposite-sign internal states, and (iii) an internal dynamics with two chaotic attractors. We 
describe the main changes in the pattern dynamics as a function of the density of elements 
in the system. A summary and a brief discusion of the results and the possible applications 
of this model are given in Sec. V. 
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II. THE SOCIAL GAS MODEL 

We propose a model of social motility that exhibits the following properties: (i) similar 
individuals attract each other; (ii) individuals can form stable groups; (iii) a group of similar 
individuals breaks into subgroups if it reaches a critical size; (iv) interaction between groups 
can modify the distribution of the elements as a result of fusion, fission, and pursuit; (v) 
individuals can change their internal state by interaction with their neighbors. We consider 
a system of motile elements where each element i is characterized at discrete time t by a 
state xl and a position vector in a (i-dimensional space. The dynamics of the system is 
given by the following set of equations. 



The evolution of both the state and the vector position of element i depends on its 
interactions with the elements in its neighborhood, rjl, such that rjl = {j : |r^ — rj| < R}. 
The parameter e expresses the coupling intensity between the states of elements, f{x) is the 
functional form that governs the internal dynamics of each element, and raj is the number 
of neighbors of the ith element. 

In Eq. ([2l), the first factor is a vector that determines the direction of motion of element 
i. Note that the influence of the neighboring elements on determining the direction of 
motion of i depends only on their angular distribution with respect to element i; the greater 
the asymmetry of the angular distribution of neighbors, the greater the magnitude of the 
displacement of i. The second factor x]^-^ J2jeri^ ^t+i gives the sense of motion, and it can be 
interpreted as the "affinity" of i for its neighborhood; if the affinity is positive, the ith element 
will move toward the denser region of the angular distribution of neighbors, otherwise, the 
element will move away from this direction. The parameter 7 expresses the coupling of 
an element with its neighborhood; if its value is too big, groups cannot be formed because 
the change ~ rj| is greater than R. In this way, the magnitude of the displacement 
depends on the position of the element relative to its vicinity and the internal states of 
its neighbors. If the element is isolated its position will be the same and its internal state 
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will be determined by f{x). This model contains a simple individual strategy based on the 
ability of individuals to evaluate their environment and react to it. In most previous models 
the influence of neighbors on a given element can be expressed as the sum of individual 
contributions Il0| . whereas in our model, due to the product of the 

two factors in Eq. ([2]), the effect of a neighbor on the motion of a given element depends 
on the states and positions of all other neighbors. This kind of holistic evaluation of the 
environmental conditions is a model property that allows the study of the social behavior of 
individuals with the ability to process the bulk of the information they acquire from their 
senses and react in consequence. 



III. CRITICAL SIZE OF A GROUP AND CRITICAL DENSITY OF GROUPS 

A stable group is deflned here as a set of elements that in isolation remain close during 
long periods of time due to the mutual interactions. If during evolution each element is a 
neighbor of the remaining elements in the group, we say that the group is a Constant Vicinity 
Group. On the other hand, if the neighborhood rjl of element i in the group is nonstationary, 
we say that the group is a Variable Vicinity Group. The diameter of a Constant Vicinity 
Group is always less than R and it can vary in time. In contrast, a Variable Vicinity Group 
pulsates in a chaotic way and the maximum distance between elements frequently exceeds 
R. 

To estimate the maximum size that each of the above types of group can reach, we 
consider a group where all the elements maintain their internal state stationary, i.e., xl_^_i = 
x\ = X. For a Constant Vicinity Group we estimate the maximum size A^^i by the condition 
\rl_^_^ — rl\ < R] that is, an element i in the periphery of the group remains a neighbor of 
the rest of elements in the group in the next time step. Assuming a uniform distribution of 
elements on a disk of diameter R, Eq. ([2]) yields 

In the case of a Variable Vicinity Group, the maximum size N2 can only be roughly 
estimated due to the non-uniform distribution of the elements. Assume that at time t the 
group is very asymmetric and the maximum distance between elements is less than R. Then, 
at time t + 1 this group splits into two or more noninteracting groups when the condition 
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7X^(A^ — 1)^ > 2R is satisfied. Therefore, tlie maximum size of a Variable Vicinity Group is 
given by 

A'. « /f -1. (4) 

In order to characterize the transition from a Constant to a Variable Vicinity Group, we 
define the quantity 

t=l i=l 

which describes the average fraction of elements that belong to a neighborhood, where 
the time average is calculated after discarding a number of transients. Figure [T] shows the 
asymptotic quantity Poo as a function of A^, with fixed values of parameters i? = 5, 7 = 0.01, 
and X = 1. For N < Ni = 27, the system forms a Constant Vicinity Group for which Poo = 1- 
At iV = A^^i the system experiences a transition to a variable vicinity group for which Poo < 1. 
At = = 31 the system undergoes a structural change, resulting in the splitting into 
two or more groups. The values of A^ at which the above transitions occur agree with the 
theoretical values of A^^i and N2 given by Eqs. [3] and HI as indicated in Fig. [H 

The temporal behavior of the system can be characterized by the distance between the 
two furthest separated elements in the group at time t, denoted by Dt. Figure [1] also shows 
the bifurcation diagram of Dt as a function of A^. For each value of A^, 5000 iterations were 
discarded as transients and the next 100 values of Dt were plotted in FiglH The size of the 
system was increased by adding one element at a random position within a distance R of the 
center of the group of size A^. We observe that Dt displays a period two for A^ < A^^i, as a 
result of the pulsatory behavior of the size of the constant vicinity group. For Ni < N < N2, 
Dt shows a chaotic behavior corresponding to a variable vicinity group. At A^ = A"2 + 1, 
the group becomes unstable and in a few iterations it splits into two or more stable groups. 
The insets in Fig. [1] show typical configurations of the system for the values of parameter 
indicated. 

Now consider the situation when A^ elements are distributed on an area of size L x L. 
When the density p = N / Li^ is below a critical value, the elements can self-organize into 
different configurations of noninteracting stable groups. In particular, configurations of 
identical stable groups of size A^^i can support the largest number of elements in the system. 
The maximum number of noninteracting groups of size A^^i that can exist in a compact 
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arrane; ement is A^^ = {L / Rf / {2^/?,) , so the mean density of elements in the system is 
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In spite of the fact that N2 determines the critical number above which an isolated group 
splits, when the interactions between groups are considered the relevant quantity to estimate 
the maximum density of elements that the system can suport in noninteracting groups is pi 
(i.e. p2 = 2(4/^ — l)/(9-\/3i?^) ~ Pi/2). However, a stable configuration with all groups 
of the same size is reached only with very particular initial conditions. For most initial 
conditions the groups in the final stable configurations have a variety of sizes. Consequently, 
in general the configurations become unstable at a critical density below pi. From now on, 
we express the surface densities of elements in the simulation area, p = N/L"^, in terms of 
the normalized density p = p/pi 

IV. SOCIAL GAS BEHAVIOR 

To illustrate the main properties of the social gas model we analyze its behavior for three 
cases of the internal dynamics f{x) in a two-dimensional system with periodic boundary 
conditions. We focus on the dependence of the pattern dynamics on the particle density p 
and show results for a particular set of model parameters, namely, R = 10 and 7 = 0.01. 

The collective behavior of the system is characterized here by the average number changes 
of neighbors per element during a number of iterations r, 

. t+T N 

i i=l 

where AnJ(+) (AnJ(— )) is the number of elements that enter (leave) the vicinity of element 
i in the time step from t to t + 1. In the following, three cases of internal dynamics f{x) 
are considered: (i) the state of all elements is stationary and homogeneous; (ii) the state 
of elements is stationary with one half of them having the state x\ = 1 and the other half 
x\ = —1; and (iii) the states of elements evolve following Eq. ([1]) with an internal dynamics 
described by the Miller- Huse map {l^: 

—2a/?i — ax, if X E (—1,-1/3) 
f{x) ={ ax, if X e (-1/3,1/3) . 

2a/3 — ax, if x E (1/3,1) 



For these three cases Fig. [2] summarizes the behavior of the system by means of the 
quantity Siooit) = Sit,r = 100) calculated each 100 iterations of t. The median Sioo,med of 
5*100 and the value 5'ioo,io% ('S'100,90%) below which S'loo fall 10% (90%) of times are shown as 
functions of p. Figure [3] shows snapshots of the system state at particular values of p. The 
main results are as follows. 

(i) The case where the state of all elements is stationary and homogeneous, i.e. f{xl) = 1, 
is shown in the top panel of Fig. [2l Note that SiQo^med exhibits a discontinuous transition 
at pc ~ 0.7 from a configuration of noninteracting stable groups with constant vicinity (as 
shown in the left top panel of Fig. [3] and characterized by S'loo = 0) to a configuration 
where elements tend to form groups with variable vicinity that become unstable due to 
sporadic interactions with neighbor groups (as shown in the right top panel of Fig. [3] and 
for which S'loo > 0). This first-order transition occurs because, as soon as one of the 
Constant Vicinity Groups exceeds the critical size A^^i, it becomes a Variable Vicinity Group 
and almost duplicates its radius. In this situation, the groups in the system can no longer 
remain isolated and they exchange elements. In general, the transition is observed at pc < 1. 
Stable configurations for p < 1 can be reached only for very specific initial conditions that 
yield to a set of isolated stable groups all having sizes < Ni. For p < Pc the system evolves 
from the disordered initial random conditions to a stable configuration like the one shown 
in the left top panel of Fig. [31 For p < Pc most initial conditions yield stable configurations 
of Constant Vicinity Groups (i.e., S = 0); however, cases with one or more isolated Variable 
Vicinity Groups can occur (i.e., S' > 0). The number of iterations necessary to reach a stable 
configuration quickly increases as p approaches pc- Reference shows 2D movies of the 
evolution of the system, corresponding to the cases shown in the top panels of Fig. [31 

(ii) The case where the state of elements is stationary with one-half of them having the 
state x\ = 1 and the other half x\ = —1 displays a behavior that is reminiscent of two 
communities that chase each other. Initially the elements are distributed at random in the 
simulation area; they soon tend to segregate into a number of groups, each group having 
members of the same sign. If the density is low enough, isolated stable groups soon form, 
but as the density increases the interaction among groups makes it difficult to reach a stable 
configuration. The existence of states of opposite sign allows for attractive or repulsive 
interactions between groups. These interactions cause the larger groups to pursue and 
disperse the smaller groups of contrary sign in frontlike configurations. At the same time, 
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attraction among elements of the same sign tends to increase the size of existing groups and 
to form new ones; however, if a group reaches a size N2 the group sphts. Groups of opposite 
signs having similar size attract each other; as a result both groups end up dispersing. As 
shown in the middle panel of Fig. S'ioo,med exhibits a discontinuous transition at 0.4. 
As stated above, for p < p^, each community forms isolated Constant Vicinity Groups 
consisting of elements with the same state (as shown in the middle left panel of Fig. [3]). 
For p > pc a chasing behavior emerges in the system (as shown in the middle right panel of 
Fig. [3]). As p increases the system displays a large degree of disorder, sustained by frequent 
interactions between elements of contrary sign that inhibit the persistence of well-defined 
groups. Movies of the evolution of the system, corresponding to the cases shown in the 
middle panels of Fig. [3] are also available at |20|. 

(iii) In the last case the states of elements evolve following Eq. ([1]) with the internal 
dynamics given by the Miller-Huse map. When a G (1,2), this map has two symmetric 
chaotic band attractors, one with values x] > and the other with x] < 0, separated by a 
finite gap centered at x = 0. For parameter values a = 1.9 and e = 0.3, the gap is small and 
interactions are sporadically strong enough to produce element sign switching. In addition 
to the chasing phenomenon observed in case (ii), now the fleeing elements may change their 
sign; a scenario reminiscent of a process of conversion or assimilation of individuals by 
larger comunities. In contrast to cases (i) and (ii), in case (iii) the isolated groups display 
randomlike motions (i.e., the chaotic evolution of the state of the elements in a group result 
in a fluctuating motion of the geometrical center of the group). For these parameter values 
the botton panel of Fig. [2]shows S'ioo,med, 'S'loo.ioro and 5*100,90% as functions of p, as well as the 
fraction fband of elements that at the final iteration {t = 6 x 10^) are in the more populated 
chaotic band, the positive or the negative. Even when the average value of the state of 
elements in a positive group is less than 1, in this case p is normalized to the same value pi 
used in cases (i) and (ii). In contrast to cases (i) and (ii), the system never reaches long- 
term stationary configurations because, even at low densities, the isolated groups move and 
eventually interact with other groups, resulting in an exchanges of neighbors (i.e., S > 0), 
and sometimes the internal state of one or more elements change their sign. As shown in the 
botton panel of Fig. [2] there is not a clear transition as in cases (i) and (ii). Note the change 
of behavior when all the elements end up in one band, i.e., fband = 1. In these simulations 
fband ~ 0.5 for 0.55 < p < 0.85. In this range of p there is exchange of elements between the 
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two chaotic bands, but an important majority is rarely established, and consequently the 
behavior is similar to case (ii). For p > 1 transient stable configurations (i.e., Siqq = 0) are 
not longer attained. The botton panels of Fig. [3] show snapshots of typical configurations 
for low and high densities (see also [20]). Figure S] shows the evolution of 6*100 ioT p = 0.4 as 
well as the fraction of elements with positive states. Note that at the beginning the number 
of positive elements is in a slight minority, but they soon become the majority, and finally 
all elements end positive. As long as fbandif) < 1, the evolution of 6*100 shows quiet periods 
(S* ~ 0) followed by active periods that are usually associated with sign changes. As soon 
as f bandit) = 1, the active periods reach higher values of S*ioo and the quiet periods are 
characterized by S*ioo = 0. This is why in the botton panel of Fig. [2] S*ioo,io% ~ when 
fband = 1 and p < 1. For p > 1 quiet periods are absent (5*ioo,io% > 0) because quasistable 
configurations are no longer possible. 

V. SUMMARY AND CONCLUSIONS 

We have analyzed the evolution of the patterns displayed by a coupled map gas in which 
the state and position of elements vary as a result of the interactions with their neighbors. 
The motion of elements is governed by a rule inspired by the fact that the reaction of a social 
individual to its environment is the result of a global evaluation instead of an evaluation of 
the sum of the neighbors' actions. The proposed model is completely deterministic, possesses 
a small number of parameters, and exhibits a series of properties that are reminiscent of 
the behavior of comunities in competition (i.e., fission, fusion, and pursuit of groups of 
elements). The pattern dynamics depend on the model parameters 7 and R and on the 
density of elements p = N/L"^ in the system area. We have analyzed the behavior of 
the system for three cases of the internal dynamics corresponding to (i) a homogeneous 
and stationary internal state, (ii) two stationary but opposite internal states and (iii) an 
internal dynamics with two chaotic attractors of opposite sign. In cases (i) and (ii) there 
is a transition at a critical value of the density from a stable configuration (quiet mode) to 
a pattern of interacting groups (confiict mode); in case (ii), pursuit and fiight of groups is 
the dominant feature in the confiict mode. In case (iii) isolated groups display ramdonlike 
displacements of their geometrical center and consequently even at low densities the system 
displays quiet periods separated by periods of confiict. The duration of the confiict periods 
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increases with p; and for high enough densities quiet periods are absent. Additionally, in 
case (iii) the sign of the elements may change and in many cases all the elements end up 
trapped in one of the chaotic bands. 

The results in this paper indicate that this kind of simple and deterministic model might 
be used to study basic properties of the collective social behavior. It must be noted that this 
study is by no means complete. Qualitatively different behaviors are expected for other sets 
of parameters and internal dynamics f{x). For example, for negative values of the parameter 
7 an isolated pair of elements repel (attract) if the elements have the same (opposite) sign; 
if all elements have the same sign stable groups cannot form, but when both signs are 
present an interesting behavior arise in which quasistable inhomogeneous groups form. If a 
periodic or quasiperiodic internal dynamic is adopted, the formation of synchronized groups 
is expected for appropriate parameter values. If the internal dynamics is multidimensional, 
as in the Axelrod model of cultural dissemination 17,[l^,[l^, the interaction among elements 



can be designed to include more sophisticated relations. 
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FIG. 1: Poo (continuous curve) and Dt/ R (dots) as functions of the group size N for parameter 
values = 5, 7 = 0.01, and x = 1. The critical sizes iVi and N2 are indicated. The insets show 
the typical configurations ior N < Ni, Ni < N < N2, and N > N2. 
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FIG. 2: Average number changes of neighbors per element as function of density when R = 10, 
7 = 0.01, and L = 50. The continuous curve corresponds to the median value -Sioo.med) whereas 
the dotted and dashed curves correspond, respectively, to the percentiles S'ioo,io% SiQo^gQ%; 
see text. Starting from random initial conditions, these quantities are calculated in the period 
10"^ < t < 6 X 10*^ and arc normalized to the critical size A^i (= 39 for these parameter values). 
Top panel: the state of all elements is stationary and homogeneous, i.e., f{xl) = 1. Middel panel: 
the state of elements is stationary with one half of them having the state = 1 and the other 
half xl = —1. Bottom panel: the internal dynamics is given by the Miller-Huse map with a = 1.9 
and e = 0.3; the dots correspond to the fraction fband of elements with the majority sign. 
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FIG. 3: Typical configurations at low density (p = 0.3, left hand panels) and intermediate density 
{p = 1, right hand-panels) at a particular iteration tg. Top, middle and bottom panels correspond 
to the three cases in Fig. [2j The positions of particles at time tg are represented with circles 
for positive states and squares for negative states. The arrows represent the displacement vectors 
between the positions at iterations — 1 and + 1. 
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FIG. 4: Evolution of Siqq normalized to the critical size A'^i (continuous curve) and the evolution 
of the fraction of elements with positive states (dashed curve). The simulation corresponds to the 
case when p = OA and the internal dynamics is given by the Miller-Huse map with parameter 
values as in the bottom panel of Fig. [2j 
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